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, Abstract. Studies of the Fundamental Plane of early-type galaxies, from small to 

^Nj ■ intermediate redshifts, are often carried out under the guiding principle that the Fun- 

damental Plane reflects the existence of an underlying mass-luminosity relation for 
^ , such galaxies, in a scenario where elliptical galaxies are homologous systems in dy- 

■ namical equilibrium. Here 1 will re-examine the issue of whether empirical evidence 

supports the view that significant systematic deviations from strict homology occur 
in the structure and dynamics of bright elliptical galaxies. In addition, I will discuss 
possible mechanisms of dynamical evolution for these systems, in the light of some clas- 
sical thermodynamical arguments and of recent N-body simulations for stellar systems 
^ . under the influence of weak collisionality. 

' 1 Introduction 

(N 

' This article focuses on three main questions: (1) Are elliptical galaxies struc- 

(«~^ I turally similar to each other? (2) Which detailed dynamical mechanisms can 

make elliptical galaxies evolve? (3) Are there general trends to be anticipated 
for the evolution of these stellar systems? 

Here I will report on a long-term research project aimed at providing answers 
O I to the above questions. Some interesting clues have been discovered only very 

^ recently |S| , , d] • Most of the paper refers to the class of bright ellipticals only; 

, low-luminosity ellipticals are known to be characterized by different dynamical 

properties. 



;_j ■ 2 Structure of bright elliptical galaxies 

, 

The answer to whether elliptical galaxies can be considered to be structurally 
similar to each other depends on the specific context in which the question is 
posed and addressed. Below, I will focus on the context of the physical interpre- 
tation of the Fundamental Plane j^Hp. 

As demonstrated by a number of investigations (e.g., see |^, (30] for a 
study based on a sample of more than 200 early- type galaxies), the observed 
correlation that defines the Fundamental Plane, logi?e — a log ctq -I- PSBg + 7 
(with a = 1.25 ± 0.1, /? = 0.32 ± 0.03, 7 = -8.895 in the B band; the effective 
radius being measured in fcpc, the central velocity dispersion in km/ sec, the 
mean surface brightness in mag/arcsec^ 0, |29jV is remarkably tight, with a 
scatter on the order of 15% in Re- 
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The following simple physical argument has been put forward as an inter- 
pretation of this important physical scaling law. If we note that (1) the ob- 
served luminosity law of bright elliptical galaxies appears to be universal (the 
so-called R^^'^ law; [221) ^^'^ (2) the kinematical structure of these systems is 
regular and uniform |2H1), it is natural to conclude that elliptical galax- 
ies should be considered as homologous dynamical systems, in the sense that 
the relevant virial coefficient Ky should be taken to be approximately constant 
from galaxy to galaxy. Then, (3) given the existence of the virial constraint, 
the Fundamental Plane can be seen as the manifestation of a mass-luminosity 
relation for galaxies (see |2ZI, IIE))- In fact, the virial theorem can be written as 
GMi,/Re = L{G / Re){Mi, / L) = Kv<Tq, where M^, is the mass of the luminous 
component and L is the total luminosity. By eliminating ctq from the Fundamen- 
tal Plane relation, one finds: 



Kv 



^ ^(2-10/3+q)/q^(5/3-q)/q ^ ^(5/3-q)/q^ 



(1) 



The latter relation follows from the empirical fact that 2 — 10/? -f a « 0. 

Unfortimately, there are empirical and theoretical findings that work against 
the hypotheses at the basis of the previous argument. First of all, significant 
deviations from the R^^^ law have long been noted (see ^Tj, 001 )j and found 
to correlate systematically with the galaxy luminosity (see also [21] )• Second, 
studies that have measured the amount and distribution of dark matter in el- 
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Fig. 1. The best-fit n(<f'), obtained by fitting tlie /oo models, projected along tlie line of 
sight, with R^^" profiles. Note the plateau at n = 4 reached by concentrated (high-!i') 
models, for which the radial range adopted in the fit is 0.1 < R/ Re < 10 (from [S]) 
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Fig. 2. The virial coefficient for the f^c (triangles) and for the isotropic R^^" (squares) 
models, based on an aperture of radius Re/8 (from [5|) 



lipticals (see HI) have shown that the presence of dark matter is more promi- 
nent in brighter and spatially larger galaxies, thus demonstrating that the virial 
coefficient may vary significantly from galaxy to galaxy. A curious theoretical 
point adds further caution to the perception that ellipticals should be consid- 
ered homologous systems. This derives from direct inspection of the so-called 
/oo sequence of models ^ . As demonstrated in |S] , models that appear to be all 
(see Fig. 1, for W in the range 7 — 10) very well fitted by the R^^'^ law, over a 
luminosity range of more than ten magnitudes, may be characterized by signifi- 
cantly different values of the relevant virial coefficient (see Fig. 2, the triangles 
representing the virial coefficient for the /oo sequence of models), as a result of 
the impact of a more and more concentrated nucleus. 

In |S] we have further confirmed, by close inspection of four cases (NGC 
1379, NGC 4458, NGC 4374, NGC 4552; studied in great detail by comparing 
the performance of a number of fitting procedures on data taken from jl9| . |18|). 
that the Sersic [Jj^ index n for the i?^/" photometric profiles can indeed be very 
different from 4 (in particular, for NGC 4552 we find n « 11, with residuals on 
the order of 0.2 magnitudes; a fit performed with the i?^/** law would lead to 
residuals up to one magnitude, while a fit based on an R^/^+ exponential proffie 
would have residuals up to half a magnitude). On the other hand, we have 
checked that, if the luminosity range where the fit to the photometric profile is 
performed is reduced to less than 5 magnitudes, then (see the profiles are 
indeed well fitted by a "universal" i?^/^ law. 
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In conclusion, while wc find it necessary to dismiss strict homology as a 
viable description of elliptical galaxies in relation to the interpretation of the 
Fundamental Plane, the existence of the empirical scaling law suggests that some 
kind of weak homology must be enforced (expressed by Eq. Q), as a correlation 
between structural properties and total luminosity. In we have also proved 
that a large scatter in the dynamical correlations (e.g., in the n ~ — 19 + SlogL 
relation noted in ^T], |^) may well be compatible with the observed tightness 
of the Fundamental Plane. 



3 Mechanisms of dynamical evolution 

Given the conclusion that elliptical galaxies have to be considered only weakly 
homologous systems, it is natural to ask whether and how individual galaxies 
may change their internal structure via dynamical processes. This general issue is 
especially important, if we recall that typically, in the study of the cosmological 
evolution of the Fundamental Plane (see and references therein), strict ho- 
mology and thus a mass-luminosity relation is assumed for the observed galaxies 
and an interpretation of the data (see Fig. 3) is made in terms of pure passive 
evolution (through the evolution of the luminosity resulting from the evolution 
of the properties of stellar populations). 

Besides the possibility of major merger events, are there significant sources 
of dynamical evolution for elliptical galaxies to be considered? As noted recently 
0, the traditional approach to the study of elliptical galaxies, in terms of equi- 
librium and stability for the solutions of the coUisionless Boltzmann equation, 
supplemented by the Poisson equation, may be misinterpreted. Given the very 
large values of typical star-star relaxation times in elliptical galaxies (see |20|. 
j42| 'l it is generally taken for granted that, unless a system happens to be in a dy- 
namically unstable state (for example, a condition of excessive radial anisotropy; 
see [33]), its state is basically "frozen" into an equilibrium distribution function. 
Thus the only task left to the dynamicist would be to decipher which distribution 
function best describes the observed states (a task that is particularly difficult 
for non-spherical systems) taken to be strictly stationary. 

In our opinion, the above picture is oversimplified and may lead to an im- 
proper perception of the dynamics of real stellar systems. If, for simplicity, we 
take the view that elliptical galaxies have formed via coUisionless collapse (see 
|45|). we should realize that splitting past and present conditions (that is forma- 
tion processes and a coUisionless equilibrium state) is just an idealization that 
the theory makes in order to define a basic state and to study its properties. In 
reality, stellar systems evolve continually and we should check to what extent 
the evolution processes change the internal structure of galaxies. 

There are several specific mechanisms and causes for dynamical evolution 
that could be studied: (i) "Granularity" in phase space left over from the ini- 
tial collapse; (ii) Presence of gas in various phases, especially of the hot X-ray 
emitting interstellar medium; (iii) Interactions with a compact central object; 
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Fig. 3. The Fundamental Plane in the rest frame B band. In panels (a) to (e), field 
E/SO galaxies are shown, binned in redshift, and compared to the Fundamental Plane 
found in the Coma Cluster by [21 • Panel (f) shows the average offset of the intercept 
of field galaxies from the local Fundamental Plane relation as a function of redshift 
(large filled pentagons) compared to the offset observed in clusters (open squares). The 
solid lines represent the evolution predicted for passively evolving stellar populations 
formed in a single burst at z = 1, 2, 5 (from top to bottom). This figure is taken from 
|44| where full references are given to the sources for cluster data points and stellar 
synthesis models 



(iv) Interactions between the galaxy and its own globular cluster system; (v) 
Interactions with external satellites and effects of tides and minor mergers. 

In a recent paper we have tried to quantify the role of items (iv) and (v) 
above by means of N-body simulations. The idea at the basis of these studies is 
that heavy objects can suffer dynamical friction and then be dragged in toward 
the galaxy center, as studied earlier, for example, in ^H], ^1], 021; in f^^ct, the 
parallel momentum transport relaxation time Tfr is related to the deflection 
relaxation time To by a factor that can be very small when a heavy test particle 
moves through a field of lighter particles: Tfr = 2TDmf/{mt + rnj). We have 
thus revisited the problem of simulating the orbital decay of a satellite, placed 
initially on a circular orbit at the periphery of a galaxy, and basically confirmed 
the general findings presented in |14j ; note that our simulations have been made 
with about one million particles, while the earlier simulations had been car- 
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Fig. 4. The development of pressure anisotropy in a galaxy as a result of the interaction 
with a shell of A'' = 100 fragments dragged in toward the galaxy center by dynamical 
friction. The broken line represents the evolving value of Kt/2, where Kt is the total 
kinetic energy associated with the star motions in the tangential directions; the solid 
line represents the evolving value of Kr, the total kinetic energy associated with the 
star motions in the radial direction (from 

ried out with a few thousand particles. Then we have proceeded to address a 
quasi-spherical problem in which the satellite is fragmented into many smaller 
objects (several runs were made with either 20 or 100 fragments), distributed 
on a spherical shell. The quasi-spherical symmetry that characterizes this study 
has the important advantage of allowing for a smoother framework, basically 
free from other effects unrelated to dynamical friction, such as those associated 
with lack of equilibrium in the initial configuration. Furthermore, with respect 
to the earlier studies of the orbital decay of a single satellite, our attention here 
is mostly shifted to measuring the evolution of the underlying structure of the 
hosting galaxies. One effect observed, while the fragments are slowly dragged in 
toward the center, is a general change in the stellar density distribution with re- 
spect to the initial polytropic basic state. Another expected effect that we have 
been able to quantify, starting from an initially isotropic distribution of stel- 
lar orbits, is the slow growth of a tangentially biased pressure anisotropy (see 
Fig. 4). All these slow dynamical evolution effects appear to be genuinely associ- 
ated with the process of dynamical friction exerted by the stars on the minority 
component of heavier objects. We are planning a survey of cases that should al- 
low us to identify general properties of dynamical evolution in elliptical galaxies 
resulting from the interaction between the stars and a significant population of 
globular clusters or of the merging of a large number of small satellites. 

4 General trends from thermodynamical arguments 



In order to study possible general trends that may be anticipated for the evo- 
lution of elliptical galaxies, we refer to the general framework that has been 
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successfully applied to the context of the evolution of globular clusters. Globular 
clusters appear to be well represented by King ^] models (see H^). They are 
recognized to be non-homologous stellar systems, subject to dynamical evolu- 
tion resulting from internal effects (such as weak coUisionality and evaporation) 
and external perturbations (such as disk-shocking, when, in our Galaxy, their 
orbits happen to cross the disk). It has been noted that these mechanisms of dy- 
namical evolution make a globular cluster evolve approximately along the King 
equilibrium sequence (see |^ and references therein). 

For globular clusters, an important paradigm is provided by the gravothermal 
catastrophe 37 ^ , which offers interesting applications and physical interpretation 
(for a review, see 02] )■ Here we recall that, starting from the study of isother- 
mal gas spheres JSl, the gravothermal catastrophe is expected to occur also in 
stellar systems (see P, [SZ])- The instability is interpreted as due to the cu- 
rious property of self-gravitating systems of being characterized by a negative 
effective specific heat. Although for stellar systems a rigorous proof has been 
provided only for idealized models where an isothermal set of stars is confined 
by a spherical box, the paradigm is generally believed to be sufficiently robust to 
be applicable to real stellar systems, provided that they possess a sufficient level 
of internal coUisionality. An independent element that strengthens the view that 
the paradigm is indeed robust has been added by an analysis that has shown, for 
an isothermal gas, that spherical symmetry is not a necessary ingredient |35| . 

Following some arguments initially put forward by Lynden-Bell (see |36| . 
|37)). would there be a way to lay out a similar scenario for elliptical galaxies 
as partially relaxed stellar systems? If so, we would gather powerful "thermo- 
dynamical" arguments to determine general trends for evolution, beyond the 
specific paths produced by a given dynamical mechanism. 

In our view, there are two aspects of the problem that require clarification. 
A first point is that we would like to start from a physically justified equilibrium 
sequence, much like King models for globular clusters, able to describe the gen- 
eral properties of elliptical galaxies. A second point is that, formally, the origin 
of the gravothermal catastrophe can be traced to the Poincare stability of lin- 
ear series of equilibria (see |22|)- For a proper mathematical derivation, one 
would thus like to start from a sequence of coUisionless models derived rigorously 
from the Boltzmann entropy. In the absence of such a sequence, a derivation of 
the gravothermal catastrophe has been based on either an unjustified ansatz (see 
1^). that the global temperature of the system would be associated with 
the coefficient multiplying the energy in the distribution function, or the use of 
non-standard entropies |U (but for unrealistic models). 

In order to address the first point, we may refer to a sequence of models that 
have been found to be very promising for a realistic description of elliptical galax- 
ies (the so-called foe models; j9j, see the review 11 ). These models have been 
inspired by the characteristics of the products of coUisionless collapse, as derived 
from N-body simulations 021 ■ In the simple spherical case, they are based on the 
distribution function f^o = A{—E)^^^exp{—aE — cJ^/2), with A,a,c positive 
constants, and define a one-parameter equilibrium sequence, which, much like 
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King models, can be parameterized in terms of the dimensionless central poten- 
tial If" — —a<P{0). For positive values of E the distribution function is taken to 
vanish. When 'F increases beyond a certain value, around 'P = 7, the models 
have a projected mass density profile that is well fitted by the R^/"^ law and 
indeed they turn out to be an excellent tool to fit the observations. From the 
point of view of statistical mechanics, they have been found 021 to be com- 
patible with a derivation based on a partition of phase space in terms of the 
star energy and the star angular momentum square, under the assumption that 
detailed conservation of the star angular momentum is required at large values 
of angular momentum. This closely follows our understanding of the process of 
partial violent relaxation . Unfortunately, the derivation is based on heuristic 
arguments and the distribution function does not follow from a straightforward 
exact mathematical extremization of the Boltzmann entropy; in particular, the 
orbit time that acts as a weight to the cells in phase space is replaced, for sim- 
plicity, by a factor l/{—E)^^^, which is approximately correct only for weakly 
bound orbits. Therefore, attempts at using this equilibrium sequence to study 
the gravothermal catastrophe in the context of elliptical galaxies, while defi- 
nitely appealing from the physical point of view (see also 8 ) , would remain less 
satisfactory from the formal point of view. 

Now we have shown |12| that we can carry out a program that is satisfactory 
not only from the physical point of view (because it is based on an equilib- 
rium sequence, also inspired by studies of coUisionless collapse [321, that is able 
to match the properties of observed galaxies), but also from the mathematical 
point of view (because the distribution is derived rigorously from the Boltzmann 
entropy by requiring the conservation of a third global quantity Q, in addition 
to total energy and total mass). The program is made possible by the second 
option explored in |43| for the construction of models of partially relaxed stellar 
systems. This option leads to the so-called /''^■' models. It was already noted 
[15] that the general physical properties of the f'^'^^ models are close to those of 
the /oo models and, in particular, that for v in the range 0.5 — 1 their projected 
mass distribution, for concentrated models, follows the R^^'^ law. 

Let / be the single-star distribution function, E the single-star specific energy, 
and J the magnitude of the single-star specific angular momentum. Consider the 
standard Boltzmann entropy: 

S = - J f\nfcPv<fx (2) 
and extremize it under the constraint that the total mass 

M ^ j fd^v(fx, (3) 

the total energy 

Etot = \j Efd?vcPx, (4) 

and a third global quantity 
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Fig. 5. Specific entropy and total energy along the equilibrium sequence of j'"^' models 
with V — 1 (as a function of the concentration parameter at constant M and Q, and 
thus expressed by means of the dimensionless functions and Note that for 

If' < 3.5 the models are characterized by a negative global temperature, because the 
derivatives of S and Etot have opposite signs. This figure is taken from \1'2] 



Q = J r\E\-^''''^fd^vd^x (5) 
are assigned. Then the resulting distribution function is 

^AeyiY,{~aE~dr\E\-'^''/^). (6) 

In the above expression, the quantities a, and d are positive constants. The 
parameter v is a, free (positive) parameter, which was argued j43| to be in the 
range 0.5 — 1.0. In the following we refer to the case v = 1. Note that the 
three constants appearing in the distribution function define two scales and one 
dimensionless parameter, which we take to be 7 = ad'^/" /{At:GA). 

Self-consistent models generated by such distribution function are computed 
from the Poisson equation, solved under the boundary conditions of regular po- 
tential at the center and of Keplerian potential at very large radii. For positive 
values of E the distribution function is taken to vanish. If we introduce the 
dimensionless central potential W = — a<?(0), the outer boundary condition de- 
fines a sort of eigenvalue problem that is solved by the relation 7 = 7(tf'). The 
self-consistent models thus make a one-parameter equilibrium sequence. 
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Fig. 6. The instability "spiral" of f^'^^ models with u — 1. The solid line refers to the 
results obtained from the ansatz that the coefficient a represents the inverse global 
temperature. Crosses represent the inverse global temperature from the definition 
dS/dEtot; other symbols indicate estimated points for which the adopted numerical 
differentiation is less reliable. Point A marks the onset of the gravothermal catastrophe 



By careful numerical integration, one may then proceed to calculate the func- 
tions S ~ S{M, Q, W) and Etot — Etot{M, Q, ^) on the equilibrium sequence (see 
Fig. 5) and from here the inverse global temperature 



The onset of the gravothermal catastrophe is thus determined by inspection of 
the equilibrium sequence studied in the [EtotX) plane (following PI]). 

In | 12| we have implemented the above program and shown that for the f^'^^ 
models the gravothermal catastrophe is expected to set in at tZ/ « 9. Surprisingly, 
around this value of the concentration, the projected mass distribution turns out 
to be very well fitted by the R^^^ law (this general point had already been noted 
in 021 ) but outside the context of the gravothermal catastrophe) . For values of 
\l/ close to and beyond 9, the general properties of the instability "spiral" in the 
(EtotjC) plane, based on the proper thermodynamical definition of the global 
temperature, are the same as in the (Etot, a) plane, based on the ansatz that the 
temperature of the models is determined by the coefficient a (see Fig. 6). 

One important point noted in jl2| is a qualitative departure of the behavior of 
the instability "spiral" at low values of 'F. For the original gas sphere and for the 
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stellar dynamical analogue of a stellar system confined by a box with reflecting 
walls, the limit of low concentration was identified as that of a non- gravitating 
ideal gas, subject to Boyle's law. In our case, the analogy breaks down. In fact, 
the global temperature turns out to change sign at w 3.5 (see Fig. 5). Such 
a drastic event should be accompanied by some physical counterpart in the dy- 
namical behavior of the system. Surprisingly, the value of If" « 3.5 coincides with 
that for the threshold of the radial-orbit instability (for the context of foo 
models, see [Z] and JD]). In other words, by undertaking a thermodynamical de- 
scription of the equilibrium sequence of models defined by the /f"^' distribution 
function, we have found arguments that lead us naturally not only to the inter- 
pretation of the observed i?^/'' law, but also to one clue for the interpretation 
of the radial-orbit instability of coUisionless stellar systems. Besides the proper- 
tics just outlined, one important additional aspect that makes the f^'^'' models, 
at this point, more appealing than the foo models is their anisotropy level. We 
had noted (e.g., see ^I]) that the foo models are actually too isotropic, when 
compared with the final products of simulations of coUisionless collapse |45| . 
The present models turn out to be much more interesting even in this respect. 
We have indeed checked that their characteristic anisotropy profile, for values 
of ^ close to the onset of the gravothermal catastrophe, is very similar to that 
observed in the numerical simulations (see Fig. 7). 




Fig. 7. Pressure anisotropy profiles a — 2 ~ {{v^) -\- {'"l))/{'Vr) as a function of radius 
for selected models {v = 1) compared to the pressure anisotropy profile found |45| 
in numerical simulations of coUisionless collapse. This figure has been prepared by M. 
Trenti 
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5 Conclusions 

For a physically justified family of equilibrium models, representing the result 
of incomplete violent relaxation, and derived rigorously from the Boltzmann 
entropy, wc have shown that, at high concentration values, the onset of the 
gravothermal catastrophe is found to occur at « 9, in the parameter domain 
where models are characterized by an i?^/^ projected density distribution. At 
low concentration values, the equilibrium sequence presents a drastic departure 
from the limit of the classical isothermal sphere, because models become as- 
sociated with a negative global temperature. The transition point, !Z' » 3.5, 
turns out to coincide with the point of the sequence where the radial-orbit in- 
stability sets in. In the intermediate concentration regime, 3.5 < < 9, the 
structural properties of the models change, much like those of models along the 
King equilibrium sequence, a family of models that is known to capture the non- 
homologous properties of globular clusters. It is our hope that, in this domain 
of intermediate concentration values, the /^''^ models may be used to describe 
the characteristics of weak homology of elliptical galaxies. 



Acknowledgements I would like to thank L. Ciotti, M. Del Principe, T. Liscikina, 
M. Lombardi, F. Pcgoraro, M. Stiavelli, M. Trenti, T. Treu, and T. van Albada, 
for their collaboration in this work, and the organizers of the workshop "Galaxies 
& Chaos" , for their warm hospitality in Athens. This work has been partially 
supported by MIUR of Italy. 

References 

1. V.A. Antonov; Vcstnik Lcningr. Univ. no. 19, 96 (1962) (Engl. Transl.: in Structure 
and Dynamics of Elliptical Galaxies, ed. by T. de Zeeuw (Reidel, Dordrecht 1986) 
pp. 531-548 

2. R. Bender, R.P. Saglia, B. Ziegler et al.: Astrophys. J. 493, 529 (1998) 

3. G. Bertin: 'Gravitational plasmas'. In Plasmas in the Universe, 142nd Course of 
the International School of Physics "Enrico Fermi", od. by B. Coppi, A. Ferrari, 
E. Sindoni (Societa Italiana di Fisica, Bologna 2000) pp. 373-393 

4. G. Bertin, F. Bertola, L.M. Buson et al.: Astron. Astrophys. 292, 381 (1994) 

5. G. Berlin, L. Ciotti, M. Del Principe: Astron. Astrophys. 386, 149 (2002) 

6. G. Bertin, T. Liseikina, F. Pegoraro: paper submitted (2002) 

7. G. Bertin, F. Pegoraro, F. Rubini, E. Vcspcrini: Astrophys. ,J. 434, 94 (1994) 

8. G. Bertin, R.P. Saglia, M. Siavelli: 'Spiraling into the R}^'^ law of ellipticals'. In 
New ideas in Astronomy, ed. by F. Bertola, J.W. Sulentic, B.F.Madore (Cambridge 
University Press, Cambridge 1988) pp. 93-95 

9. G. Bertin, M. Stiavelli: Astron. Astrophys. 137, 26 (1984) 

10. G. Bertin, M. Stiavelli: Astrophys. J. 338, 723 (1989) 

11. G. Bertin, M. Stiavelli: Rep. Prog. Phys. 56, 493 (1993) 

12. G. Bertin, M. Trenti: Astrophys. J. 584, 729 (2003) 

13. W.B. Bonnor: Mon. Not. Roy Astron. Soc. 116, 351 (1956) 



Weak Homology of Bright Elliptical Galaxies 



13 



14. Tj.R. Bontekoe: Orbital decay of satellite galaxies. PhD Thesis, Groningen Uni- 
versity, Groningen (1988) 

15. Tj.R. Bontekoe, T.S. van Albada: Mon. Not. Roy. Astron. Soc. 224, 349 (1987) 

16. A. Burkcrt: Astrou. Astrophys. 278, 23 (1993) 

17. N. Caon, M. CapaccioU, M. D'Onofrio: Mon. Not. Roy. Astron. Soc. 265, 1013 
(1993) 

18. N. Caon, M. Capaccioli, M. D'Onofrio: Astron. Astrophys. Suppl. 106, 199 (1994) 

19. N. Caon, M. Capaccioh, R. Rampazzo: Astron. Astrophys. Suppl. 86, 429 (1990) 

20. S. Chandrasekhar: Astrophys. J. 97, 251 (1943) 

21. P.H. Chavanis: Astron. Astrophys. 386, 732 (2002) 

22. G. De Vaucouleurs: Ann. d'Astrophys. 11, 247 (1948) 

23. S. Djorgovski, M. Davis: Astrophys. J. 313, 59 (1987) 

24. S. Djorgovski, G. Meylan: Astron. J. 108, 1292 (1994) 

25. M. D'Onofrio, M. Capaccioli, N. Caon: Mon. Not. Roy. Astron. Soc. 271, 523 
(1994) 

26. A. Dressier, D. Lynden-Bell, D. Burstein et al.: Astrophys. J. 313, 42 (1987) 

27. S.M. Fabcr, A. Dressier, R.L. Davies et al.: 'Global scaling relations for elliptical 
galaxies and implications for formation'. In Nearly normal galaxies: From the 
Planck time to the present, ed. by S.M. Faber (Springer, New York 1987) pp. 
175-183 

28. O. Gerhard, A. Kronawitter, R.P. Saglia, R. Bender: Astron. J. 121, 1936 (2001) 

29. I. J0rgensen, M. Franx, P. Kjairgaard: Astrophys. ,J. 411, 34 (1993) 

30. I. J0rgensen, M. Franx, P. Kjaergaard: Mon. Not. Roy. Astron. Soc. 280, 167 
(1996) 

31. J. Katz: Mon. Not. Roy. Astron. Soc. 183, 765 (1978) 

32. J. Katz: Mon. Not. Roy. Astron. Soc. 189, 817 (1979) 

33. J. Katz: Mon. Not. Roy. Astron. Soc. 190, 497 (1980) 

34. I.R. King: Astron. J. 71, 64 (1966) 

35. M. Lombardi, G. Bertin: Astron. Astrophys. 375, 1091 (2001) 

36. D. Lynden-Bell: Mon. Not. Roy. Astron. Soc. 136, 101 (1967) 

37. D. Lynden-Bell, R. Wood: Mon. Not. Roy Astron. Soc. 138, 495 (1968) 

38. M. Magliocchetti, G. Pucacco, E. Vesperini: Mon. Not. Roy. Astron. Soc. 301, 25 
(1998) 

39. V.L. Polyachenko, I.G. Shukhman: Sov. Astron. 25, 533 (1981) 

40. P. Prugniel, F. Simien: Astron. Astrophys. 321, 111 (1997) 

41. J.L. Sersic: Atlas de galaxias australes (Observatorio Astronomico, Cordoba 1968) 

42. L. Spitzer: Dynamical evolution of globular clusters (Princeton University Press, 
Princeton 1987) 

43. M. StiaveUi, G. Bertin: Mon. Not. Roy. Astron. Soc. 229, 61 (1987) 

44. T. Treu, M. StiaveUi, S. Casertano et al.: Astrophys. J. Lett. 564, 13 (2002) 

45. T.S. van Albada: Mon. Not. Roy Astron. Soc. 201, 939 (1982) 

46. T.S. van Albada, G. Bertin, M. StiaveUi: Mon. Not. Roy. Astron. Soc. 276, 1255 
(1995) 

47. E. Vesperini: Mon. Not. Roy. Astron. Soc. 287, 915 (1997) 

48. M.D. Weinberg: Mon. Not. Roy. Astron. Soc. 239, 549 (1989) 



